Finite nuclear size and Lamb shift of p-wave atomic states 
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We consider corrections to the Lamb shift of p-wave atomic states due to the finite nuclear size 
(FNS). In other words, these are radiative corrections to the atomic isotop shift related to FNS. It is 
shown that the structure of the corrections is qualitatively different from that for s-wave states. The 
perturbation theory expansion for the relative correction for a j>i/2-state starts from a\n(l/Za)- 
term, while for Si/2-states it starts from Za 2 term. Here a is the fine structure constant and Z 
is the nuclear charge. In the present work we calculate the a-terms for 2p-states, the result for 
2pi/2-state reads (8a/97r)[ln(l/(Za) 2 ) + 0.710]. Even more interesting are p 3 / 2 -states. In this case 
the "correction" is by several orders of magnitude larger than the "leading" FNS shift. 

PACS numbers: 11.30.Er, 31.30.Jv, 32.80.Ys 



Experimental and theoretical investigation of the ra- 
diative shift (Lamb shift) of energy levels in heavy atoms 
is an important way to test Quantum Electrodynamics 
in presence of a strong external electric field. One of 
the effects related to this problem is a dependence of the 
Lamb shift on the finite nuclear size (FNS) . One can also 
look at this effect from another point of view. It is well 
known that there is an isotop shift of atomic levels due 
to the FNS. The corrections we are talking about are the 
radiative corrections to the isotop shift. 

The corrections for Is-, 2s-, and 2p-states have been 
calculated numerically, exactly in Za, in Refs.|l], |^, [|. 
The self-energy and the vertex corrections to the FNS 
effect for any s-wave state have been calculated analyti- 
cally in order a(Za) in Refs.|Q, [|. However, the struc- 
ture of the higher order Za corrections and, in partic- 
ular, their logarithmic dependence on the nuclear size 
has not been understood even for s-states. Our interest 
to FNS radiative corrections has been stimulated by our 
work on the radiative corrections to atomic parity non- 
concervation . Technically the parity nonconservation 
effect has some common features with that of the FNS 
radiative correction: in both cases the effective size of 
the perturbation source is much smaller than the Comp- 
ton wavelength Ac- In the paper || we have elucidated 
the structure of higher order in Za FNS radiative cor- 
rections for s-electrons, and have calculated analytically 
a(Za) and a(Za) 2 ln(Ac/Vo) self-energy and vertex FNS 
relative radiative corrections. Here ro is the nuclear ra- 
dius. In the present work we calculate FNS radiative 
corrections for p-wave electrons. We demonstrate that 
the structure of the corrections for p-wave states is very 
much different from that for the s-wave states. Physically 
it happens because of different infrared behavior. 



Due to the finite nuclear size, the electric potential 
V(r) of the nucleus is different from that for a pointlike 
nucleus. The deviation is 



SV(r) = V(r) - 




Throughout the paper we set h = c = 1. The dia- 
gram that describe the FNS effect in the leading order 
is shown in Fig.|l](a). The double line corresponds to the 
exact electron wave function in the Coulomb field, and 
the zigzag line with cross denotes the perturbation (nh. 
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FIG. 1: The leading contribution to the FNS effect is given by 
diagram (a), and one loop radiative corrections to the effect 
are given by diagrams (b-e). The double line denotes the exact 
electron Green's function in the Coulomb field of the nucleus, 
the cross denotes the nucleus, the zigzag line denotes the FNS 
perturbation ^), and the dashed line denotes the photon. 

Diagrams Fig.|l](b) and Fig.^(c) correspond to the con- 
tributions of the electron self-energy operator and the 
vertex operator, respectively. The diagram Fig.|l](d) de- 
scribes a modification of 6V (see eq. (Q)) due to the vac- 
uum polarization, and the diagram Fig.[l](e) corresponds 
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to a modification of the electron wave function due to the 
polarization of the vacuum by the Coulomb field (Uehling 
potential). 

Technically the most complicated are the self-energy 
and the vertex FNS (SEVFNS) corrections given by di- 
agrams in Fig.[l](b) and Fig.|l|(c). According to our pre- 
vious work H, the SEVFNS relative correction for an 
s-wave state is of the form 
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Here b = exp(l/(2 7 ) - C - 5/6), C w 0.577 is the Euler 
constant, and ro, as we already mentioned, is the nu- 
clear radius. The total relative SEVFNS correction (||) 
is the ratio of the sum of diagrams FigJ^(b) and Fig.|l|(c) 
divided by the diagram Fig.|l|(a). Value of A s is not pro- 
portional to the nuclear radius squared because it is a 
relative quantity. Plot of A s versus the nuclear charge Z 
is shown in Fig.|| by the dashed line. Results of computa- 
tions of A s for Is and 2s states Q are shown by squares 
and triangles, respectively. The a(Za) term in (0) comes 



is no contribution of larger distances into A s ? The rea- 
son is very simple. In the leading order the correction A s 
can be expressed in terms of the forward electron-nucleus 
scattering amplitude H. There is a rigorous QED theo- 
rem that claims that there is no an infrared divergence 
in the forward scattering amplitude, see e.g. Ref. [||. 
Therefore, quantum fluctuations from distances r Ac 
cannot contribute to A s (see also Ref. ||). Let us look 
now at the p-wave SEVFNS correction A p . From the 
point of view of the scattering problem it corresponds 
to scattering at finite angle. The finite-angle scattering 
amplitude is always infrared divergent. Therefore, one 
must expect a contribution to A p from quantum fluctu- 
ations at distances r 3> Ac- This is the contribution we 
calculate in the present work. 

Formally we assume that Za <C 1. Therefore, at dis- 
tances r ~ Xc/(Za) 3> Ac dynamics of the electron is de- 
scribed by usual nonrelativistic Coulomb wave functions. 
However, the nucleus radius is small, ro <C ZaXc -C Ac- 
At so small distances, generally speaking, one must use 
relativistic Dirac wave function even at Z = 1 . The elec- 
tron Dirac wave function at r <C ZaXc is of the form 



tf(r) = Nr~<- 
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FIG. 2: Relative SEVFNS corrections (%) for s 1 / 2 - andpi/ 2 - 
states. The dashed line shows the correction A„ given by En. 
(0). Results of computations of A s for Is and 2s states W 
are shown by squares and triangles, respectively. The solid 
line shows the correction A p given by Eq. Jj"?|). Results of 
computations of A p forlp^^ state j^/ are shown by diamonds. 

from distances r ~ Ac, and the a(Za) 2 ln(Ac/?"o) term 
comes from distances ro <C r <C Ac- An important point 
is that there is no contribution that comes from distances 
r ^> Ac- Because of this reason the correction A s is ex- 
actly the same for Is, 2s, 3s,... states |7|. Why there 



where VL and = — (er • n)fl are spherical spinors ; k — 
— 1 for si/2-state, k — 1 for p!/ 2 -state, and k = —2 for 
p 3 / 2 -state ; 7 = \J ' k 2 — (Za) 2 ; and N is a constant 
known for each particular state, see Ref ||. For s 1 / 2 - 
and p 3 / 2 -states the upper component of the Dirac spinor 
(^) is much larger than the lower one. Hence, the upper 
component determines the FNS shift of such a state. On 
the other hand, for p]/ 2 -state the lower component and 
hence its contribution to the FNS shift is dominating. A 
straightforward calculation gives the following values for 
the FNS shifts of 2s and 2p states (diagram FigQa)) 



— (Za) 4 m 3 <r 2 >, 



SE (0 \2 Pl/2 ) = ^(Zafm 3 < r 2 > , 
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(5£ (0) (2p 3/2 ) = -L(Za) 6 m 5 < r 4 > 
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Here < r 2 > and < r 4 > are values of r 2 and r 4 averaged 
over charge density of the nucleus. The low- momentum 
expansion of the nuclear electric form factor is of the form 



F(g2)s , 1 _^ <r 2 >+ _L_ <r 4 > (5) 

Modeling the nucleus as a uniformly charged ball one gets 
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where ro = 1.1 A l / 3 fm is the nucleus radius, and A 
is the nucleus mass number. As one should expect 
the FNS corrections (0) obey the following inequalities 
5E(°)(2s 1/2 ) » 6EM(2p 1/2 ) » 5EM(2p 3/2 ). 

Let us calculate now the leading in Za one loop 
SEVFNS radiative correction E^ for pi/ 2 - and P3/2- 
states. This correction is given by diagrams in Fig.|l](b) 
and Fig.|l|(c). Since we consider the leading correction, it 
is sufficient to use the nonrelativistic approximation for 
electron wave functions (two-component wave functions) . 
It is sufficient also to use the effective FNS perturbation 
that reproduces FNS correction for s-wave states, 



Sy eff(r) = gS(r) 
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Rest of the calculation is very similar to the textbook 
calculation of the Lamb shift, see, e.g. Ref. ||. We in- 
troduce the parameter \i such that m 3> (i m(Za) 2 . 
Hence the correction E^ 1 ' can be represented as a sum 
of "high frequency" and "low frequency" contributions 
E^ = Eh + Ei, where "high" and "low" correspond 
to frequencies above and below (X, respectively. In the 
momentum representation, the effective potential corre- 
sponding to the high frequency contribution is of the form 
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where q = Pi — p 2 is momentum transfer, and 7 is the 
Dirac matrix. Taking the p-wave component of the po- 
tential (^|) and transferring it to the coordinate represen- 
tation, we get the following expression for the SEVFNS 
high frequency correction for a p-wave state \ip) 
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This gives the following values for 2p-states 



E h (2p 1/2 ) = F 



E h (2p 3/2 ) 
where F 
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The contribution of the vacuum polarization, diagram 
Fig.|l|(d), can be taken into account in Eqs. (p|),(|l0|) by 
substitution ln(m/2/z) — > \n{m/2(i) — 1/5. The Uchling 
potential, diagram Fig.|l|(e), does not contribute in this 
order. 

The low frequency contribution is given by the usual 



nonrelativistic quantum mechanics expression 
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Here 5i stays for real part, lj is frequency of the virtual 
photon, H = p 2 /2m — Za/r is the nonrelativistic Hamil- 
tonian, and E^ = —m(Za) 2 /& is the energy of 2p-state. 
We have also taken into account that interaction with the 
photon is of the form —ep ■ A/m, where A is the vector 
potential of the photon. The contribution Ei is the same 
for pi/ 2 - and for p 3 / 2 -state. Using explicit form of 2p 
wave function, one can represent (jllj) as 
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where 
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Eigenvalues of H are e n = —m(Za) 2 /2n 2 . Therefore, 
the first impression is that the integrand in Eq. ([l2]) is 
singular at lj = and u) = e 2 — e%. However, the function 
4>{r) is orthogonal to the wave function il> 2s {r), hence, 
there is no real singularity at u = 0. There is a real 
singularity at u> — e 2 — ei that is related to the possibility 
of emission of real photons, and this slightly complicates 
integration in (|l2]). To overcome this technical problem, 
it is convenient to represent <j>(f) as (f>{r) = ip(r)+/3ipi s (r) 
with (3 = y / 2~77r(2/3) 4 . In this form <p(r) is orthogonal 
both to ip 2s {r) and tpi s (r). Then ( p^ ) is transformed to 
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where |0) 



In 



denotes the electron localized at origin, 
this form the matrix element (0| (e^ — H — ui\ \<p) 



has no singularities. Using explicit expression for the 
nonrelativistic Coulomb Green's function llOl 



G(0,r\E) = (r\ 
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where p = \/—2mE, r\ = mZa/p, T is the gamma- 
function and W is the Whittaker function, and taking 
the integral over r, and then over u>, we obtain 



Ei=F 
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Combining ([To]) and (|l6|), we finally obtain the total 
SEVFNS radiative corrections (diagrams Fig.[j](b) and 
Fig.llVc)) in the leading order 



E^(2 Pl/2 ) = F 
E^(2p 3/2 ) = F 
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As one should expect, the result is independent of the 
parameter fi. We have already mentioned that to ac- 
count for the vacuum polarization (the diagram Fig.[l](d)) 
one has to replace ln(l/(Za) 2 ) ln(l/(Za) 2 ) - 1/5. 
Therefore the total FNS radiative corrections (diagrams 
Fig.[l](b), Fig.[l](c), and Fig.[j](d)) in the leading order are 



E£l(2 Pl/2 ) = F 
E^ t (2p 3/2 ) = F 
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Let us have a look now at the relative FNS radiative 
correction for 2p 3 / 2 state. According to Eqs. (||) and 
ll8|) the relative correction is 

20 a < r 2 > 
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For example, for Hydrogen atom the radiative correction 
is by a factor 2.6 10 4 larger than the "leading" contribu- 
tion. 

According to the present calculation, the leading in 
powers of Za SEVFNS relative radiative correction (di- 
agrams Figjl](b) and Fig.|l](c)) for 2p!/ 2 -state is equal to 
(see comment [|Ll| ) 
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As we have already explained, the correction comes from 
quantum fluctuations at distances A c < r < X c /(Za). 
There is also a contribution 



a(Zaf 
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that comes from distances r <C r <C \c, this contribu- 
tion has been calculated in our previous work || . The 



contribution cx a(Za) that comes from r ~ Ac has not 
been calculated yet. Therefore, altogether one gets the 
following formula for the relative correction A p : 



- a 



- — In 

9vr I (Za) 2 



0.910 



(21) 



ax(Za) 



(Za) 



where a\ is an unknown coefficient. To determine the 
coefficient a\ , we fit results of numerical calculation of A p 
for 2p 1 / 2 -state [Q. As a result of the fit we find a± = 2.75. 
The correction A p given by Eq. (^l|) is plotted in Fig.^ 
by the solid line. The results of computations 0] are 
shown by diamonds. Agreement is very good. 

Concluding, we have shown that corrections to the 
Lamb shift of p-wave atomic states due to the finite nu- 
clear size are qualitatively different from that for s-wave 
states. The difference is related to the infrared behavior 
of quantum fluctuations. As a result, the leading relative 
p-wave correction is proportional to a hi(l/Za) while the 
leading s-wave correction is proportional to a(Za). The 
leading p-wave correction has been calculated analyti- 
cally. 
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